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We study the oscillations of neutron and strange stars in gravity. More precisely the nonradial 
/-modes are examined and the differences with pure general relativity are investigated. Using these 
results we build several gravitahonal wave asteroseismology relahons. Our goal is to determine up 
to what extend these relations are equation of state independent and whether they deviate enough 
from general relativity in order to produce an observable effect. The results show that the differences 
coming from R^ gravity are up to 10% and that will be difficult to be observed in the near future. On 
the other hand the small deviahons in some of the asteroseismology relations show that they are not 
only equation of state independent, but they are also quite insensitive to the gravitational theory. That 
is why solving the inverse problem can give us quite robust estimates of the neutron star parameters. 

PACS numbers: 


I. INTRODUCTION 


In the last few years the accelerated expansion of the universe has been confirmed, leading fo the investigations of 
different cosmological models capable of explaining it. A new type of matter should be introduced when considering 
this problem in the context of general relativity (GR) - the so-called dark energy which constitutes about 73% of the 
total energy content of fhe Universe and exhibifs some exotic properties, like a negative pressure-to-density ratio. 
A way of overcoming fhis is fo explore the generalized theories of gravity as a possible alternative explanation for 
fhe accelerafed expansion. One of the most popular representatives are the f(R) theories which allow us to exclude 
the dark energy hypothesis. These theories of gravity are a natural generalization of GR, derived by replacing the 
Lagrangian in the Einstein-Hilbert action with a more general one lll-lsl]. 

The viability of any theory of gravity should be tested against the astrophysical observations too. In this paper 
we are investigating one such astrophysical application for a specific type of /(R) theory, namely the so-called 
R—squared gravity, where the form of the Lagrangian we adopt is /(R) = R + aR^. Differences between general 
relativity and alternative theories are expected to occur for strong gravitational fields, such as the ones created by 
different compact objects like neutron stars (NS), strange stars (SS) and black holes (BH). The scope of our studies is 
limited to the first two of them. 

Taking in consideration that the new generation of gravitational wave detectors will start operating in the next few 
years, we extend the studies of compact stars in /(R) gravity to gravitational wave asteroseismology. The neutron 
star perturbations can lead to nonradial oscillations, which are a source of gravitational waves. It is expected that 
in the near feature the emitted gravitational waves will be observed by the earth based detectors and the obtained 
information about the oscillation frequencies will reveal a new spectrum of research possibilities for astrophysics. 
The full development of gravifational wave asteroseimology will let us use optimally the gained information. The 
relations and the methodology that would allow us to obtain the neutron star parameters via the gravitational wave 
observations have been developed in the last two decades in many papers for G R M -lgl and for some alternafive 
models iiS, and in the last few years even for rapidly rotating neutron stars lfl^.Eolk 
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More detailed comments on some of the above-mentioned papers, which our work is based on, should be made. 
In Is , 211 the authors investigated the fimdamental and the pressure mode oscillation frequencies and damping times 
as functions of the average density \/M/R^ of the stars or as functions of the compactness M/R for equations of 
state (EOS) with different stiffness. They also introduced some empirical relations between these quantities that 
can be used to infer the neutron star mass and the radius, i.e. solve the inverse problem. Later those relations 
were reexamined for a set of newer EOS by Benhar et al. Lau et al. @] suggested a different normalization, 
namely to use the so-called “effective compactness" rj = VM^/ 1, with 1 being the moment of inertia, instead of the 
compactness M/R, motivated by the studies in ll23l . Using that parameter they introduced an empirical relation 
valid for both NS and SS, which turned out to be quite EOS independent. In a recent paper Chirenti et al. @] 
investigated oscillations of neutron stars in GR with all of the above mentioned relations. They commented on 
the EOS independence of these empirical relations. Other EOS independent relations cormecting the neutron star 
moment of inertia, quadrupole moment and the tidal Love numbers were considered for example in I23l- i3^ . 

For the first time the oscillation frequencies of neutron stars in alternative theories of gravity, and more specifically 
in the case of a scalar-tensor theory (STT), were investigated in BEi. There one can also find a detailed derivation 
of the perturbation equations for STT in the so-called Cowling approximation (the background metric and the scalar 
field are fixed). 

The relations presented in the above-mentioned papers are EOS independent to a large extent. Such relations are 
a useful tool for overcoming the uncertainty naturally caused by the EOS and can be used to explore generalized 
theories of gravity. This will be our main goal in the present paper, namely to explore the various asteroseismology 
relations in R?- gravity and to quantify the deviations from general relativity. Such results can be potentially used 
when solving the inverse problem, i.e. obtaining the stellar parameters from the observed oscillation frequencies, in 
order to impose constraints on the gravitational theory. In the current paper we will concentrate on the nonradial 
oscillations for which the pressure is the restoring force. The scope of our study is limited to the fundamental modes 
(f-modes) that are the major object of investigation in the literature cited above. The reason is that the pressure 
modes (p-modes) are a much less efficient source of gravitational waves, so it is expected that the main contribution 
to observational data will be from f-mode oscillations. 

It is important to point out that in some of the above-mentioned papers, and more precisely in the more com¬ 
plicated cases concerning rapid rotation and alternative theories of gravity, the Cowling approximation is used 


1 15l. Il9l. 120(1 . It is known that there are differences between the results obtained via the Cowling approximation 
and the full general relativistic results, and the deviation decreases with the increase of the compactness. However, 


19t|20|]. That is why we will employ it in our 


this approximation is accurate enough for qualitative investigations 
studies. 

The structure of our paper is as follows: In Section II the reduced field equations describing neutron stars in R?- 
gravity are presented, in Section III the equations for the perturbations are examined, in Section IV the results are 
presented and discussed. The paper ends with conclusions. 


II. BASIC EQUATIONS 

In this section we briefly present the basic equations describing equilibrium neutron and strange star solutions in 
R?- gravity. More details on this problem can be found in i3ll,l33l . 

The f{R) theories are described by the following action: 

S ~ I67tG J 5matter{g}iv,X)i (1) 

where R is the scalar curvature with respect to the spacetime metric gfiv The action for the matter fields, denoted 
by Xr is separated in the term Scatter- The following inequalities have to be satisfied for the f{R) theory to be free of 
ghosts and tachyonic instabilities: 


dR^ 


> 0 , 


dR 


> 0 . 


( 2 ) 
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In the case of gravity, which we are considering, the above inequalities give a > 0 and 1 + 2aR > 0. 

A useful fact, that we will employ in our studies, is the mathematical equivalence between f{R) theories and 
Brans-Dicke theory (with cobd = 0 and nonzero scalar field potenfial) given by the action 

S = J d^x^g [OR - 11(0)] + Smamr{g^v,x), (3) 

where O and the potential L1(0) are defined as follow: O = and 11(0) = R^ — f{R). In fhe case of R^ gravity 
0 = 1 + 2aR and the Brans-Dicke potential is L1(0) = ^ (O — 1)^. 

While the above action is in the so-called Jordan, or physical, frame, for mathematical simplicity it is useful to study 
the scalar-tensor theories in the so-called Einstein frame. The metric in this frame is defined by the conformal 
transformation = Og^i, and the action can be written in the form 

s = ^ J [R* - 2g*f^^d^cf>dyf - V{cp)] + SMatter{e~^^^glyX), (4) 

where R* is the scalar curvature with respect to the Einstein frame mefric g’^y. The new scalar field is defined by 

(p = ^ InO, the potential in the Einstein frame is given by V[q>) = and A^{cp) = 0”^((p). The field 

equations in this frame are much simpler and fhat is why we will use it. Since the physical quantities are measured 
in the Jordan frame, we will make fhe fransition between the two frames where necessary. 

In our sfudies we will consider oscillafions of stafic neutron stars. But in some of the asteroseismology relations 
the moment of inertia is used. A natural and straightforward way to define this quantity is via the slow rotation 
approximation. Eor this reason we will briefly examine here the more general framework of slowly rotafing neutron 
stars in R^ gravity. 

The line element in a stationary and axisymmetric spacetime, keeping only first-order terms in the angular velocity 


Q, can be written in the form: 

dsl = + r^{dd^ + sin^ Qdd^) — 2(v{r, 6)r^sin^Qdddt. (5) 

The explicit form of fhe field equafions in the Einstein frame is as follows 

^ 87rGA4(<p)p + + ^v{cp), (6) 

(^) -^V{cp), (7) 

$ - (S - ^ ^ - 4.C.,- 3.,- . 

f = -(p+p)(^+«w5). m 

^ A ^ = i(,jiGA‘^{cp){p +p)a>{r), (10) 

r* dr I dr j 

where we have defined 

^(^) ^ dlm^(^ and (d = Cl-cv. (11) 


Here p and p are fhe energy and pressure density in the Jordan frame and they are connected to the Einstein frame 
ones via p* = A^{(p)p and p* = A^(^)p. 

The above sysfem of equations, combined with the equation of stafe for the matter and appropriate boundary 
conditions, describes the interior and the exterior of a compacf sfar. The boundary conditions are the natural ones. 
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At the center of the star we have p{0) = pc,A{0) = 0, ^(0) = 0, ^(0) = 0, where pc is a free parameter denoting 
the central value of the density. The condition for ^(0) = 0 ensures the regularity of fhe mefric function lo at the 
center of fhe star, the condition ^ (0) =0 ensures the regularity of the scalar field q), and in turn the regularity of the 
Jordan frame scalar O af fhe center of the star. The regularity of the metric functions at r = 0 requires A(0) = 0 and 
since the Einstein and the Jordan frame metrics are conformally equivalent, this condition ensures also the regularity 
of the Jordan frame geomefry af r = 0. The boundary conditions at infinity are related to the fact that we consider 
asymptotically flat space-time. At infinity we have limr->oo <J’(r) = 0, limr-^ood) = Cl and lim^-^oo <J>(r) = 0 with 
V (0) = 0. These conditions ensure the asymptotic flatness for both frames. 

Although the above equations are in the Einstein frame, the final results we present in this work are in the Jordan 
frame (the physical one). The coordinate radius of the star is determined by pjrg) =0 while the physical radius 
in the Jordan frame is given by Rg = A[(p(rg)]rg. 

The explicit form of the conformal factor A{(p) and the potential V{(p) for R^ gravity are 

A^{<p) = e-T3‘P, « = = (12) 

The moment of inertia I of the compact star is defined in the standard way 



(13) 


where J is the angular momentum of fhe start. More convenient for numerical calculations is to use an integral 
expression for the moment of inertia. 


I = 


8nG 


f ^ A^{(p){p + p)e 
Jo 




(14) 


Erom now on we shall use fhe dimensionless parameter a —>• a/Rq and the dimensionless moment of inertia 
I —> I/MqRq where Mq is the solar mass and Rg is one half the solar gravitational radius Rg = 1.47664 km. 


III. NEUTRON AND STRANGE STAR OSCILLATIONS IN COWLING APPROXIMATION 


In this section we will present the equations governing the non-radial oscillations of spherically symmefric com¬ 
pact stars in the Cowling approximation for R^ gravity. In this case we are investigating fluid perfurbafions on a 
fixed Jordan frame mefric which is equivalent to fixed scalar field and mefric in fhe Einstein frame. Despite these 
simplifications the approximation turns out to be accurate enough for qualifative investigation. The results differ 
from GR in the range of 10 fo 30% and fhe deviation decreases wifh the increase of fhe compacfness ]3^ 341. 

The Jordan frame equations describing the perturbations in the Cowling formalism are obfained by varying the 
equation for fhe conservation of fhe energy-momenfum fensor in the Jordan frame, namely V^STy = 0. The Jordan 
frame Lagrangian fluid displacement vector ^ can be parameterize by two functions W and V in the standard form 




e-^W.-vC,- 


1 / 3 


sm' 


2 0 3d 


^icot 




(15) 


where the Jordan frame metric function and the Jordan frame radial coordinate r are given by = A(^)e^ and 
r = A{q))r. Equivalently ^ can be written in the form 


^ = A-^{cp){A{cp)e-'^W,-V-,- 


V 3 
0 3d 


sm' 


^icot 


^y,^(0,d). 


(16) 







5 


Then the system of ordinary differential equations describing the star oscillations in Cowling approximation is the 
following: 


d]N 

dr 

dV 

dr 


^ {a}^A[(p)e^-^‘^Vr^ + O'W + ci{cp)ipW^ - l{l + l)A{(p)e^V, 


dp 

2 (o' + a.{cp)xp) - 


l\ e^WA-'^icp) 

^ o / 


(17) 


where we have defined ip = dq)/dr. 

The boundary condition at the star surface is that the Lagrangian perturbation of the pressure vanishes, which is 
equivalent to 


+ (O' + ix{(p)xp) e-^W 


-Aw^ (<?) 


= 0 


(18) 


r=R 


From a numerical point of view it is convenient to introduce new variables Wi and Vj and express the boundary 
condition at the center of the star in terms of them 


W = W-ir^+\ V = Vir‘. 


(19) 


We are rewriting the system (flTt and the boundary condition at the surface in therms of Wi and Vi. The new 
system of equations has the following form: 


dWi 


dr 

dVi 

dr 


dp 
dp . 


(xrA{cp)e^ ^Vir + + 0i{(p)ipWi 


l{l + l)A{cp)e^Vi 


( 20 ) 


2 (o' + a{cp)ip) - 


Vi- 


e^WiA-\(p) 


with the boundary condition at the surface of the star 


+ (O' + a{cp)ip) = 0 


r=R 


( 21 ) 


and at the center 


W, = -lA{cp)V^l^,. ( 22 ) 

The above system l l20ll combined with the boundary conditions forms an eigenvalue problem with the frequencies 
CO being the eigenvalues of the system. In the next section we present the results from the numerical solution of the 
eigenvalue problem. 


IV. RESULTS 

We investigate numerically the fundamental mode (f-mode) frequencies obtained from the system of eqs. (l20ll 
supplied with six realistic hadronic equations of state and two quark ones. The eigenvalue problem is solved using 
a shooting method where the oscillation frequency plays the role of the shooting parameter. The two classes of 
equations of state cover a very wide range of masses and radii as shown in Fig. [T] EOS FPS and BBB2 are the softest 
ones with masses below the observational limit of 2M© Nevertheless we use them in order to check the EOS 
independence of the relations we derive ^. On the other hand EOS SLy, APR4 and WFF2 have masses above 2MO 


^ It should be mentioned that inf (R) theory of gravity some soft EOS could be reconciled with the observations because the maximum mass 
increases, as demonstrated in Isill . 
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FIG. 1: The mass-of-radius relations for all EOS in pure GR. It could be seen that they cover a wide range of stiffnesses. 




FIG. 2: Results for one hadronic EOS, namely SLy, and one quark EOS, namly SQSB60. Both have maximum mass around two 
solar masses. In the left panel the normalized frequency as a function of the neutron star mass is presented while in the right 
panel - the oscillation frequency in kHz as a fimction of the average density of the star in km^^ . 


and fall into the preferred range of radii 13711 . EOS MSI is the stiffest one with maximum mass over 2.7M© and much 
bigger radii. 

For the hadronic EOS we use piecewise polytropic approximations 13811 . The quark ones on the other hand have 
the analytical form 


p = b{p- po), (23) 

where the constants b and po ^re taken from for EOS SQS B60 and SQS B40. The first one gives masses slightly 
below two solar masses and the second one is stiffer and has higher masses and radii. 

We will examine a wide range of relations between the /-mode oscillation frequencies and the stellar parameters 
such as mass, radius and moment of inertia. The particular choice is motivated by the most commonly used relations 
in the literature on gravitational wave asteroseismology. First we will investigate the neutron and strange stars 
oscillation frequencies and the differences between GR and gravity. We thoroughly investigated the f-modes for 
all EOS, but here we present results for one hadronic and one strange matter EOS, namely EOS SLy and SQS B60, 
that can be considered as representative examples. The sequences of models range from one solar mass up to the 
maximum one in accordance with the observations. In Fig. |2]we present and compare the results for GR and 
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FIG. 3: In the left panel we present the frequency as a function of the average density for all examined EOS for the case of GR and 
for gravity with parameter a = 10^, which induces a deviation close to the maximum possible. It could be seen that the results 
for the stiffest EOS, namely MSI and SQSB40, seriously differ from all the other results. However, the two groups of results form 
two narrow bands. In the right panel we plot the relation Mcv{M/R). The presented results are again for all EOS for GR and 
/(R) theory with a = 10^. As one can see, the results for hadronic and quark EOS are only shifted from each other, but there is 
quite good universality for both cases independently. The separation between the GR and /(R) dependencies is negligible. 


gravity for neutron and strange stars. In the left panel the normalized frequency (jj{R ?/as a function of the 
mass of the star, measured in solar masses, is presented for EOS SLy and SBSB60. Such a relation was used in other 
studies of alternative neutron star models (iB - iilii. For clarity, the SQSB60 results are marked with additional 
symbols. The GR case is presented by a continuous black line and different values of a are in different colors and 
patterns. It can be seen that for both GR and /(R) gravity the results for neutron stars are much steeper than those 
for strange stars and the normalized frequency decreases with the increase of the mass. The maximum deviation 
from GR is below 10% and the behavior is qualitatively different for the different mass ranges. For smaller masses 
the deviation increases when a < 10 but it starts to decrease for larger a. For bigger masses though the deviation 
increases monotonically with the increase of a. 

In the right panel of Fig. |2]we present the dimensional frequency / = w/ {27z), measured in kHz as a function 
of the average density of the star (M/in km“^. This is one of the most standard asteroseismology relations 
used in the literature, originating from j^- A qualitative difference between the results for neutron stars in GR and 
in /(R) theory can be seen on the graph. The dependences in /(R) theories are no longer linear for large values of 
a. For strange stars the graphs have similar shapes for both theories. For small values of a the mode frequencies in 
R^ gravity are always higher than GR, and a significant decrease with respect to GR is observed for larger average 
densities and large values of the parameter a. 

In the left panel of Fig. |3] a graph similar to the one in the right panel of Fig. |3 is presented but in this one all 
EOS are included for GR and R^ gravity with a = 10^. This value of a is chosen because our investigations show 
that it gives almost the maximum possible deviation from the pure Einstein's theory. Different EOS are marked with 
different symbols, the GR case is in black and a = 10^ is in color. Two groups can be distinguished in the graph - one 
for the softer and the typical EOS (with maximum mass around and below 2M0), and one for stiffer EOS (contains 
EOS MSI and SQS B40). We will concentrate on the former group which could be separated into two bands, one 
formed by the results for GR and the other for a = 10^. These two bands fully overlap for small values of the 
average density, and partially for higher densities. We should note that strange stars do not fit very well in the band 
formed by neutron star solutions in Fig. iBil. Let us comment on the deviations coming from the modification of 
GR. Although the difference between GR and R^ gravity is non-negligible, i.e. up to 10%, it is comparable with the 
uncertainty of the EOS. Constraining further the EOS by different observations will reduce the spread of the data 
due to the EOS uncertainty. This will make the difference between the two theories of gravity clearer and that can be 
potentially used to set constraints on the /(R) theories using future gravitational wave observations. 
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FIG. 4: The relation Mco(t]) is presented in the graphs. In the right panel we investigate the deviations in this relation for EOS 
SLy and several different values of the parameter a. In the left panel results for all EOS are presented.The relation for both the 
neutron and the strange stars are quite EOS independent and a visible difference between the GR and the a = 10^ cases exists. 
The maximum observed deviation is close to 10%. 


The empirical relations describing these dependences have been widely investigated in the literature starting from 
|@]. They have used linear fitting for the frequency-average density relation and we will employ the same. We have 
excluded from the relation the very stiff nuclear EOS MSI and the quark ones SQSB60 and SQSB40 since they clearly 
lead to significantly different dependences and they are not favored by observations. The results for GR and for 
gravity with a = 10 ^ are fitted separately using a linear fit of the form 

[m 

/ = Ci + C2^^. (24) 

The dimensions of the constants for this fit, Ci and C 2 , are as follows: Ci is in kHz and C 2 is in kHz/km. For the case 
of GR we have C\ = 1.59, C 2 = 24.23, and for R? gravity with a = 10^ we have C\ = 1.95, C 2 = 14.25. From these 
numbers it is obvious that the results in GR are much steeper than these in / (R). 

Let us proceed to the next asteroseismology relations we plan to consider, where different normalization of the 
quantities is used. In the right panel of Fig. |3]we present the scaled frequency Moo as a function of the compactness 
of the star M/R, as proposed in j^. The graph contains all EOS, the GR solutions are marked with black squares 
and the ones for a = 10^ with circles. As one can see normalizing the frequent seriously decreases the spread of 
the data due to the EOS and therefore it leads to significant EOS independence ISlIbl kol] . The difference between GR 
and R^ gravity is below 5% which is not large enough to be observable. 

For small values of M/R the relation is almost linear, but close to the maximum values of the compactness this 
changes. That is why we use a cubical polynomial fit in this case. 


Ma; = Ci + C2- + C3(^-j • (25) 

We excluded both quark EOS and fitted the rest of our data for three different cases. For the pure GR case Ci = 
4.95 X 10-3, C 2 = 4.16 X 10-1, C 3 = 5.17, C 4 = -14.43, and for the R^ gravity Q = -2.54 x 10-3, C 2 = 6.48 x 
10-1, C 3 = 3.67, C 4 = —11.57. If one uses the data for both theories we have Cj = 2.76 x 10-^, C 2 = 0.544, C 3 = 
4.37, C 4 = -12.95. 

As one can see the results for strange stars are also quite EOS independent but they differ from the neutron star 
ones. More precisely they are shifted with some constant value of the normalized frequency. Lau et al. |@] suggested 
that this is due to the difference in the density profiles of strange and neutron stars. That is why they used a new 
parameter for constructing the universal asteroseismology relations. Instead of the compactness they employ the 
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parameter t] = where I is the moment of inertia of the star. Therefore we have also investigated the relation 

Mco{rj). In the left panel of Fig|4]we present results for EOS SLy for GR and a few different value of the parameter 
a. The biggest deviation is around 10% and it is relatively the same for all values of a > 10. In the right panel we 
present our results for all EOS for GR and for a = 10^. As we pointed out in Eig. |3l if we use the compactness M/R 
as a dimensionless parameter and normalize co, the deviation between the two theories is imder 5%. Using rj as a 
parameter preserves the EOS independence but the results for GR and gravity get further apart. However the 
deviation between GR and R^ gravity is still small, up to roughly 10% as we commented. In this case the dependences 
for both theories are more or less shifted with some constant value of Mco. As far as the difference between neutron 
and strange stars is concerned, it decreases a lot if we use rj as a dimensionless parameter instead of M/R, but the 
results in the two cases do not fully coincide as opposed to the ones in |@]. This is most probably due to the fact that 
we are not solving the full nonlinear Einstein equations, but the Cowling approximation is employed instead. 

In order to obtain asteroseismology dependences for these normalizations, we are using a quadratic fit similar to 

0 , 


Mlo = Cl + C2^ + (26) 

Eor GR we obtain Ci = —3.28 x 10“^, C 2 = 0.426, C 3 = 0.157 and for R? gravity with a = 10^ we have Ci = 
-4.29 X 10-2, C 2 = 0.544, C 3 = -2.57 x 10-2. 


V. CONCLUSIONS 

In this paper we investigated neutron and strange stars oscillations in GR and R?- gravity. More precisely we 
concentrated on the f-mode oscillations since they are expected to be very efficient gravitational wave emitters. 
We employed a big range of hadronic and strange matter EOS with different stiffness and the calculations were 
performed in the Cowling approximation. The observed maximum deviation between the f-mode frequencies in GR 
and R 2 gravity is up to 10 % and depends on the value of the R 2 gravity parameter a. 

We investigated multiple gravitational wave asteroseismology relations available in the literature and obtained 
the corresponding analytical fits. In most cases the dependences in GR and R2 gravity are qualitatively the same - 
they are only shifted with respect to one another. The only exception is the relation connecting the f-mode oscillation 
frequencies to the average density of the star where some qualitative differences between the two theories exist. As 
mentioned above these differences do not exceed 10%. Such deviations are big enough to make a clear distinguish- 
ment between GR and r 2 gravity, but as far as the real gravitational wave observations are concerned, they will most 
probably be below the accuracy of the observed gravitational wave frequencies. 

However this result also can be considered as a strong point of the asteroseismology relations examined in the 
present paper. The reason is that one can make the conclusion that the relations we consider turned out to be not 
only EOS independent, but up to a large extent theory independent too. Therefore solving the inverse problem will 
supply us with unique values for the parameters of a star, like mass and radius, insensitive to both the EOS and the 
particular theory of gravity. In order to prove this conjecture more rigorously one of course has to examine these 
relations is other alternative theories of gravity, but the present results give a hint in this direction. 
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